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Few-body correlations in the QCD phase diagram 
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Abstract. From the viewpoint of statistical physics, nuclear matter is a strongly correlated many-particle 
system. Several regimes of the QCD phase diagram should exhibit strong correlations. Here I focus on 
three- and four-body correlations that might be important in the phase diagram. 

PACS. 21.65.+f Nuclear matter - 21.45. +V Few-body systems 
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1 Introduction 

Lattice and effective model calculations provide a rich and 
exciting sketch of the phase diagram of Quantum Chro- 
modynamics (QCD). Two regions may coarsely be dis- 
tinguished: a hadronic phase and a plasma phase. Since 
quarks are considered part of the fundamental building 
blocks of matter the existence of a hadronic phase is al- 
ready an indication of strong correlations between quark- 
antiquark and three-quarks. Further on, in both phases 

, correlations lead to more interesting phenomena, such as 
clustering of nucleons to form nuclei, or superfluidity (in 
nuclear matter) and color superconductivity (in quark mat- 
ter) . In particular in the later cases a weak residual inter- 
action is enough to destabilize the ground state (just as is 

. the case for the formation of Cooper pairing). These in- 
vestigations are usually based on the study of two-particle 

\ correlations. There are reasons to go beyond two-particle 
correlations, e.g.: 

— Particle production even in a dense environment such 
as deuteron formation in a heavy ion reaction, need a 
third particle to conserve energy- momentum H . 

— To study the properties of a-particlcs (2) or determine 
the critical temperature of a possible a-particle con- 
densate |^.^|5| needs an in- medium four-body equa- 
tion. 

— Recent results in the Hubbard model indicate, that 
three-particle contributions may lead to a different (lo- 
wer) critical temperature compared to the simple Thou- 
less criterion Question of this type have not been 
addressed for nuclear matter. 

— The chiral phase transition is often discussed along 
with a confincmcnt-dcconfincmcnt transition based on 
investigating mesons (quark-antiquark states), see e.g. 
Ref. |Q . Does this transition happen for nucleons (three- 
quark states) at the same density /temperature? 

To investigate these issues a first step is to develop and 
solve proper effective in-medium three- and four-body equa- 
tions that are valid at finite temperatures and densities 



analogous to the Feynman-Galitskii or Bethe-Goldstone 
equations ||. This has been achieved in the past for the 
nonrelativistic problem (|,||,|lo],[ll]]. These equations have 
been derived on the basis of statistical Green functions Q . 
The Green functions have been decoupled utilizing a clus- 
ter expansion, see e.g. To tackle these questions in 
the (deconfincd) quark-phase, in addition, such in-medium 
few-body equations must obey special relativity. As chiral 
symmetry breaking is presumably restored (up to small 
current masses) the quarks may become very light ob- 
jects and therefore relativistic effects should play a larger 
role than for isolated systems. Here, relativity is realized 
using the light front form of relativistic dynamics [fl3|| . 
First results are given in Refs. on the confinemcnt- 

deconfinement (Mott) transition. 



2 Theory 

We use Dyson equations to tackle the many-particle prob- 
lem, see e.g. Ref. This enables us to decouple the hi- 
erarchy of Green functions. The Dyson equation approach 
used here is based on two ingredients: i) all particles of a 
cluster are taken at the same global time ii) the ensemble 
averaging for a cluster is done for an uncorrelated medium. 
The resulting decoupled Green functions may be econom- 
ically written as resolvents in the n-body space, where 
n = 2, 3, 4, ... is the number of particles in the considered 
cluster. 

The solution of the one-particle problem in Hartree- 
Fock approximation leads to the following quasi-particle 
energy 



A. 

2mi 



+ ^V 2 (12,12)/ 2 



2mf 



+ r HF (0). (1) 



The last equation introduces the effective mass that is a 
valid concept for the rather low densities considered here 
and /j, oS = (i— I7 HF (0). The Fermi function = f(ei) for 
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the i-th particle is given by 

/(e< 



1 



e /3( £i -^) _|_ I " w 

The resolvent Go for n noninteracting quasiparticlcs is 

n 

Go(z) = (z-H )- 1 N = R (z)N, = (3) 

where Go, H , and iV are formally matrices in n particle 
space. The Matsubara frequency z\ has been analytically 
continued into the complex plane, z\ — > z ^]. The Pauli- 
blocking for n-particles is 



A r = /"l/ 2 .../n±/l/ 2 .../„, 



1=1- f (4) 



where the upper sign is for Fermi-type and the lower for 
Bose type clusters. The full resolvent G(z) is given by 



G(z) = (z-H - Vy'N, 



V 



Note that V' / V. For the two-body case as well as for a 
two-body subsystem embedded in the n-body cluster the 
standard definition of the t matrix leads to the Feynman- 
Galitskii equation for finite temperature and densities H , 



T?{z) = V 2 a + V 2 a N?R (z)T?(z). 



(6) 



Introducing the Alt Grassberger Sandhas (AGS) |L6| tran- 
sition operator U a p(z) the effective inhomogeneous in- 
medium AGS equation reads 

U a(3 (z) = (1 - 5 aP )R^\z) + £ N]T^(z)R Q (z)U^(z). 

(7) 

The homogeneous in-medium AGS equation uses the form 
factors defined by 



to calculate the bound state ipB 3 

\F a ) =J2Sc t pN 2 3 Tg(B 3 )R (B 3 )\F f3 ). 

Finally, the four-body bound state is described by 
1^) = y j KrUlJB i )Ro{B i )N^{B i )R {B i )\TZ 



(8) 



(9) 



T7 



where a C a, 7 C r and a, r denote the four-body parti- 
tions. The two-body input is given in (^|) and the three- 
body input by (Jfy. Note that, although we have man- 
aged to rewrite the above equations in a way close to the 
ones for the isolated case, they contain all the relevant in- 
medium corrections in a systematic way, i.e. correct Pauli- 
blocking and self energy corrections. The numerical solu- 
tion requires some mild approximations that are however 
well understood in the context of the isolated few-body 
problem. 



3 Results 

An experiment to explore the equation of state of nuclear 
matter is heavy ion collisions at various energies. Here 
we focus on intermediate to low scattering energies and 
compare results to a recent experiment 129 Xe+ 119 Sn at 
50 MeV/A by the INDRA collaboration JTJ|. A micro- 
scopic approach to tackle the heavy ion collision is given 
by the Boltzmann equation for different particle distri- 
butions and solved via a Boltzmann Uehling Uhlenbeck 
(BUU) simulation Jl8,19|. The reaction rates appearing 
in the collision integrals are a priori medium dependent. 
However, previously this medium dependence has been ne- 
glected. Within linear response theory for infinite nuclear 
matter the use of in-medium rates leads to faster time 
scales for the deuteron life time and the chemical relax- 
ation time as has been shown in detail in Refs. pfj|.^Tt. 
This faster time scales should have consequences for the 
freeze out of fragments. 

We use the in-medium AGS equations (0) that repro- 
duce the experimental data in the limit of an isolated 
three-body system. For details on the specific interaction 
model see Ref. B . We investigate the influence of medium 
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Fig. 1. BUU simulation of the deuteron formation during the 
central collision of 129 Xe+ 119 Sn at 50 MeV/A. 
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Fig. 2. Ratio of proton to deuteron numbers as a function 
of cm. energy. The experimental data are from the INDRA 
collaboration. 
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Temperature T= 10 Me V, cm. momentum P =0 
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Fig. 3. Difference between the pole energy of the bound state 
and the continuum, B(n,T) = -E po i c — E cont - 



dependent rates in the BUU simulation of the heavy ion 
collision as compared to use of isolated (i.e. experimental) 
rates. Figure [I] shows that the net effect (gain-minus-loss) 
of deuteron production becomes larger for the use of in- 
medium rates (solid) compared to using the isolated rates 
(dashed). The change is significant, however, a compar- 
ison with experimental data is difficult since deuterons 
may also be evaporating from larger clusters that has not 
been taken into account in the present calculation so far. 
The ratio of protons to deuterons may be better suited 
for a comparison to experiments that is shown in Figure 
[|. The use of in- medium rates (solid) lead to a shape closer 
to the experimental data (dots) than the use of isolated 
rates (dashed). In these calculation, besides the change of 
rates, also the Mott effect has been taken into account. 

Figure || shows the dependence of the binding energy 
for different clusters at a given temperature of T = 10 
MeV and at rest in the medium. 

In Figure ^ part of the phase diagram of nuclear matter 
is shown. The condition for the onset of superfluidity for 
ct-particles is B(T C , fi, P = 0) = 4jU. The critical tempera- 
ture found by solving the homogeneous AGS equation for 
fi < confirms the onset of a condensation even at higher 
values (dotted) as given earlier (solid, from ||) which was 
based on a variational calculation using the 2+2 compo- 
nent of the a particle. For /i > the condition E = 4/i 
for the phase transition can also be fulfilled. However, the 
significance for a possible quartetting needs further inves- 
tigation. Due to the many channel situation of more than 
two-particles equations, the Thouless criterion might be 
revisited. 

Work supported by Deutsche Forschungsgemeinschaft. 
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Fig. 4. Critical temperatures of condensation/pairing leading 
to superfluid nuclear matter. For an explanation see text. 
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